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The Elliptic Cylinder Function of Class K. 

By William H. Butts. 



The object of this paper is the synthetic treatment of the Elliptic Function 
of Class K and the computation of tables of values that may be useful to physicists 
in studying the properties of elliptical membranes and elliptical cylinders. The 
literature on this subject is very limited, as the only serious attempts to discuss 
this function have been made by Mathieu and Heine. The investigation of 
Mathieu in his Memoire sur le mouvement vibratoire d'une membrane de forme 
elliptique [Jour, de IAouville) made some progress in the theoretical treatment 
of the function, but did not lead to definite conclusions. Heinrich Weber, in 
Annal. von Clebsch und Neumann, Bd. I, says of Mathieu's work: "Die Inte- 
gration ist dort durch Reihen bewerkstelligt, von denen mit grossem Fleisse eine 
betrachtliche Anzahl Glieder berechnet sind, fur welche aber ebenfalls kein all- 
gemeines Gesetz angegeben ist. Diese Untersuchungen mogen daher fur den 
Physiker immerhin von grossem Werte sein, mathematisch scheint mir das 
Problem dadurch der Losung wenig naher gebracht zu sein, als durch die Auf- 
stellung der gewohnlichen Differentialgleichung selbst." 

E. Heine, in Kugelfunhtionen, Bd. I, II, makes greater progress in the 

analytic treatment, but does not give a satisfactory proof* of the convergence of 

the series and does not carry the investigation far enough to make the results 

useful to the physicist. 

3 2 V d z V dPY 
Laplace's equation ~— a + ~— - % -f- ^~s =0 is reduced by Heine, after various 

transformations and the introduction of Lame's elliptical coordinates, to the 
convenient form 

* Since the above statement was written, a satisfactory proof has appeared in the Inaugural Dissertation 
of Simon Dannacher, Zurich, 1906. 
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Our investigation will be limited to the function of the Class K, 

1 °° 

(2) E ==■ ^ a + 2 a n cos 2n$, 

in which a is independent of <£> but a function of z and of the argument b. Under 
what conditions is this a solution of (1)? 

Substituting this value of E in (l), we find the following recursion formulas, 
showing the relation between the coefficients : 

(3) a x = — £ bza , 

a 2 = b (1 — z)a x — a , 
a 3 = b (4 — z) « 3 — a 1 , 



a«+i = b(n z — z) a n 



To determine the necessary conditions for the convergence of (2), it is 
necessary first to show that ^i 1 ^ a n = 0. 

Substituting a finite portion of (2) in (l) gives 

(4) j^ + (| cos 2$ + 4z S ) E = | (a n cos (2n + 2) $ — a^ cos 2n$) . 

If values of z can be found that will make a n and a n _ x arbitrarily small, 
the first member may be made to approximate zero. 

Assuming that a = 1, a ly a 3 , . . . . , a„ are rational integral functions of z, 
and equations (3), in the form a n = f(z) = 0, can be solved for various values 
of the argument b. Equations (3) give the following values of a lf a s , . . . . , a 7 , 
explicitly in terms of b and z. Placing these equal to zero we compute the 
values of z that will make a x , . . ■ . , a 7 approximately zero. 

(5) a, = — \bz = 0, 

a s=-l h *(* - 5z2 + ( 4 - p) Z + p) = °> 
a4=+ ^(^_ 1 4 2 3 + (49-i)^ + (-36 + f), + (-^ + ^)) = 0, 




a =0 

n 
ft-te. — . r 

b=0./ 
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« 5 = - \V(/- ^ + (273 _ *,) z 8 +(_ 820 + ™ 5 ) z 3 

. /_„„ 652 . 5\ , /1152 40\\ 
+ ^576-- p - + p ) 2 + (-^— v ))=0, 

a s =+ lb*(tf - 55z 5 + (l023 - |) a* + (- 7645 + 2 * 3 4 ) z 3 

+ (21076 - ™» + |) * + (- 14400 + 1™ - £) . 
, / 28800 . 1072 2\\ 

a,=z—±V (f - 91s 6 + (3003 - £) z 5 + (- 44473 + ^) ^ 

+ ( 39 62 9 6- l2 ^ + ^)^ + (--773136 + l 3 ^ 4 -^) 

+ (578400-^ + ^-l 6 ) Z + (l^-^ + ^)) = 0. 

In order to study these equations and to deduce from them properties of the 
equation a^ = 0, we compute tables of values showing the roots to four decimal 
places and the values of the functions a lf a 3 , ...., a 7 for these approximate 
roots. The Tables I, II, III are given for b = 0.1, b = 10 and b = 100. The 
last figure of every result is given exactly as found, and is not increased by unity 
if the next figure is five or more. 

Tables I, II, III and Figs. 1, 2 show how rapidly these roots approach a 
constant limiting value as « increases without limit. This line of argument 
forces us to use a^ as a function of z of infinite degree and to treat a x =0 as 
an equation of infinite degree. This extension of the function concept may be 
justified by the same necessity which forces us in certain problems to use infinity 
as a limit. 

From equations (5) we deduce the following : 



b = 0.1. 



(6 ) ^ = _ _0ji ( 8 ) = o, 



a 2 = + iML 2 (z 3 — z — 200) = 0, 

« 3 = — iML (z 3 — 5z 3 — 296z + 800) = 0, 
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a i= z +l2J! 4 ( z 4_14z 3 — 351z 3 + 3564z+ 12800) =0, 

a B = — ^^ (z 6 — 30z 4 — 227z 3 + 9680z 3 — 146 24z — 284800) = 0, 



« 6 = -f (°- 1 )' i (z 6 _ 55z 5 -f 423z 4 + 16755z 3 — 221524z 2 — 275600z 



a 7 = — QH. (z 7 — 91z 6 + 2303z 5 + 4527z 4 — 802004Z 8 + 6731264z 2 



+ 5840000) = 0, 
6731264Z 2 
+ 17224000z — 181760000) = 0. 



b = 10. 



(7) a 1 = -H (z) = 0, 



« 2 = + i^-(z 2 -z-0.02) = 0, 

in 3 

a 3 = — if- (z 3 — 5z 2 + 3.97z + 0.08) = 0, 

a 4 = +1^ (z 4 — 14z 3 + 48.96Z 3 — 35. 64z— 0.7198) = 0, 

a B = _ 1 ^ B ( 2 b_30z 4 + 272. 95z 3 — 818 .95z 2 4- 569.4805z 4- 11 .516) = 0, 

a 6 = + — (z 6 — 55z 5 + 1022. 94z 4 — 7642. 56z 3 -f 21042. 7409z 2 

— 14225. 1401z— 287.892802) = 0, 

a T = — — (z v — 91z 6 4- 3002. 93z B — 44468. lz 4 + 296172. 1714z 3 
z 

— 771755. 623z 3 + 511811 .455993z 4- 10364.025712) = 0. 



b = 100. 



(8) a 1 = - l f(z) = > 



a % = + —(z 3 - z - . 0002) = 0, 

a 3 =— — (z 3 — 5z 2 + 3.9997z+ 0.0008) = 0, 

a 4 = +l^(z 4 — 1433+ 48.9996Z 2 — 35. 9964z— 0.00719998) = 0, 

Z 



10 10 , 



a B = _i™(z 5 — 30z 4 + 272.9995Z 3 — 819.9895z 2 + 575.93480005Z 

4- 0.1151996) =0, 
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in 13 
a 6 = + ±T-(z 6 — 55z B + 1022.9994s 4 — 7644.9756a 3 + 21075. 66740009z 2 

A 

— 14398. 25120201Z— 2.879989280002) = 0, 
a 7 = — -£- (z 7 — 91z 6 + 3002. 9993z B — 44472. 951z 4 + 296294. 76170014z 3 

— 773122. 1956063z 2 + 518334. 105689599993z 
+ 103.679602560112) = 0. 

These equations (6), (7), (8) are solved by Horner's method, and the last 
remainders are multiplied by the factors before the parentheses in the equations 
under consideration. By z m we represent the third root of a 6 = 0, and by a & the 
value of a 6 for this approximate z 63 . See Tables I, II, III. 

Maximum and Minimum Values OFa w . 

The maxima and minima values of a n can not be found by the usual 
methods when n increases beyond all limits ; since a w can not be expressed 
explicitly in terms of z, and a! ^ = is an equation of an infinite degree. 

Substituting a n = b n (n — 1) ! 2 /?„,* we compute the values of z in (3' n = 
for finite values of n, z ni being taken graphically as the abscissa corresponding 
to the maximum or minimum value of (3 n between the i th and i + 1 th roots of 
(3 n — 0. We compute 5 from /?'„ = by Horner's method, accurate to one 
decimal place. Substituting these values of z in (3 n by Horner's method, we 
obtain close approximations for (3 n . 

I. Argument 6= 10. 

( 9 ) P n== b n (n—1)!*' 

(10) fa = + l(z*-z-0.02), 

/? 3 = — ^ (z 3 — 5s 2 + 3.97z + 0.08), 

& = +y2(2 4 — 14z 3 + 48.96z 3 — 35. 64z — 0.7198), 

/? 6 =— — ^(z 5 - 30z 4 + 272.95z 3 - 818.95z 2 + 569.4805z i • 11.516), 

*See Dannacher, p. 16. 



Butts : The Elliptic Cylinder Function of Glass K. \ 35 

&= + 28800 ( z6 ~ 55z5 + 1022.94z 4 — 7642. 56s 3 + 21042.7409s 3 

— 14225.1401Z— 287.892802), 

^~~ 10^600^"" 91z6 + 3002.932 5 — 44468. lz 4 + 296172.1714Z 3 

— 771755.623z 3 -f 511811.455993z + 10364.025712). 

From these equations Table V is computed. 

Substituting (9) in the recursion formula, 

(11) a n = b((n — lf~ z) a„_! — a n _ 2 , 

b n (n — 1) / 3 (3 n = b ((n — l) 2 — 2) 6"- 1 (n — 2) / 3 — 6 K - 2 (n — 3)/ 3 (3 n _ 2; 

W P» = {} — (n — iy) P*- 1 ~ #»(n— l) 8 (n— 2) 8 ^ n ~ 2 " 

The first minima for /?„ are found by computation to have the relation 
I Asi I > I |^ 3 i I >• I (3 a I , and all are negative. See Table V. 

(13) ft, = (l - ~) {3 Z1 - 10 , x 1 fl i x j r Ai, for * - + 0.4. 

(14) & 1= (l -^f)&i- 108X 1 3 8 X2 » &*' for *« = + 0.4. 

(15) & = (l - ^|) j3 fl - ^^.^ &i, for i„ = + 0.4. 

(16) ?, 1 = - %0 ^00 -x, , -lo^^T 3 P- --• i» f « 5 • 1 ■ 

The first two roots of $„ = must lie between — 2 and 0, +1 and +4, 
respectively, as proved by Heine.* Hence the first minimum lies between 
— 2 and +4. Beginning with /3 B1 , the parenthesis of (15) cannot vary from 
unity by more than one-fourth, and the parenthesis rapidly approaches the limit 
unity in ft nl as n increases without limit. The last term in (15) equals 
— 0.0007 /3 31 and in subsequent equations rapidly vanishes, so that the first term 
in (15) and in subsequent equations controls the sign and (3 nl approaches the 
limit @ n _ h 1 for z nl . Hence | /3»i | < | /?„_!, 1 1 • • • • < | fti 1 and 

(!') n=™ Pnl = Pm-l,l. 

From computations it is evident that the first minimum of fi^ = is 
numerically less than 0.11 and negative ; approximately, — 0.09. 



* Kugelfunktionen, I, 412. 
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By similar reasoning, the first maximum (5^% is +0.33, approximately. 
By further computation and reasoning from the recursion formula, it is found 
that the succeeding minima and maxima of @ n converge to limits for w=oo, 
each larger than the preceding numerically, and all differing from zero. 

The general form of the curve {3^ =0 is shown by Pig. 3, the nullpoints 
being determined by z w lf z^ 2 , . . . • , z^ 7 , as subsequently computed and given 
in Table II. 

II. Argument b = 100. 

For argument b = 100, Z?^ = has nullpoints still nearer 0, 1, 4, 9, . . . . , 
(n — l) 2 , as shown in Table III. The locus of fi^ =0, for 5=100, has the same 
general form as Fig. 3, the first minimum and the first maximum being nearly 
the same as for 5 = 10. In no case is a maximum or minimum zero. 

III. Argument b = 0.1. 

The maxima and minima values of (3^ for 6 = 0.1 present greater diffi- 
culties, due to the fact that (0.1) 3 occurs in all values of {3 n and to the fact that 
the roots of (3 n = do not fall in the regular intervals — 2, 1, 4, 9, ...., 

(n — l) 3 , re until 

(18) b(n — 2)>1. 

See Table I and Heine's Kugelfunhtionen, I, 407. 

To determine the laws governing these maxima and minima values, we 
compute z and /? found in Table IV, using the following equations : 

(19) & = + \ (s 2 - z - * ) = + I ( s 2 - z - 2 ) , 

(20) & = — I (s 3 — 5s 2 — 296s + 800), 

8 

(21) &= + £- (s 4 — 14s 3 — 351s 2 + 3564z + 12800), 

(22) (3 & = — — !-- (s B — 30s 4 — 227s 3 + 9680s 2 — 14624s — 284800), 

1152 

(23) /? 6 = + .„' (s 8 — 55s B + 423s 4 + 16755s 3 — 221524s 2 — 275600s 

28800 

+ 5840000), 

(24) ft = — ~-4f~-(s 7 — 9 Is 6 + 2303s 5 + 4527s 4 — 80 2004z 3 + 6731264s 2 

1036800 

+ 17224000s— 181760000), 
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(25) #= + l( z _i) = 0, 

(26) / 8i = _|( s ?_3i Z -98f)=0, 

(27) $ = + ± (z 3 - 10| z 3 - \1b\z + 891) = 0, 

(28) PL=— ™- (z 4 — 24z 3 — 136iz 2 + 3872z — 2924*) = 0, 

1152 

(29) #= + — *— (z»— 45fz 4 + 28 2z 3 -f 8377|z 3 — 73841£z — 45983*) = 0, 

(30) $ = — - * ann (z 6 — 78z 5 + 1645z* + 2301|z 3 — 3437162 2 + 1923218fz 

+ 2460571$) = 0. 

In Fig. 4, we locate the nullpoints of (3^ =0 from the values ofa aol toa aol 
in Table I, and draw the curves representing /3 2 = to /3 7 = 0. 

In studying these curves the following tendencies should be considered (see 
Table I and Fig. 4) : 

a) The nullpoints of these curves always approach the nullpoints of (3^ = 
as n increases. 

b) The first z nl and the corresponding minima values of (3 always increase 
numerically with n. „ L i™ (3 nl must be examined. 

c) The first maximum increases from /3 3 = to /? 4 = and afterwards 
decreases, apparently to a small positive limit for (3^ = 0. 

d) The other minima and maxima between two successive roots of /?„ = 
always decrease numerically as n increases, apparently toward a very small 
limit for {3^ = 0. 

e) In each curve, the maxima and minima for the first half of the arches 
retrograde and for the last half advance beyond the middle of the interval. 

To discuss these tendencies and eventually to discover properties of (3^ = 0, 
return to the formula 



a n = 6 ((n — l) 3 — a) a n _ x 



a 



n — 2, • 



Substituting a n = b n (n — 1) / 3 /3 n gives 
(31) (3 n = (l— (n _^ 1)8 ) ft,-i -^_ 1)3(n „ 2 ) 3 ^- 2 - 
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First Minimum. 

To explain tendency b) and to determine the limit of the first minimum, 
use equation (31) and Table IV. The first minimum evidently lies between 
z aol = — 16.9015 and z^^^ — 5.0524, the first two roots ofa cc = and of 

Beginning with /3 6 = 0, 

(32) ft, = (l + %i) p a - 0t01 *# x $ Ai , for i* = - 12.6 , 

(33) £ 61 = ( 1 + %J-) n - - 0m01 xP X 4 > P« » f0r ** = ~ 13 -° ' 

(34) ? n = (l+ I | i ) At - o.o^Vxb' ^ ' fOT ** = - 13 ' 3 ' 

(35) ? sl = (l + {5^) A-u - .01(n-l)»(n-2)» ^~ 2 ' f0r "■' ' 

Since /? B1 , /8 61 , /? n are known by computation to be negative and 

I&1KI&1KI&1I, it is evident from (35) that |&il< I AxK • • • •< !£«, i|, 
and that they are all negative, since the first term in the second member of (35) 
has the greater multiplier from /3 B1 to fi^ lf the last multiplier being a proper 

fraction and decreasing rapidly while the first multiplier 1 + , nl L g remains 
greater than unity. 

The relation | /3 n | < | /? 81 1 . . . . < | (3^ x | is preserved when z^ x is substituted 
for i n , %, . . . ., and the inequalities are still greater, since less ordinates are 
substituted for maximum ordinates. 

(36) However, 

iaki ( j + j^if)p^\ and ift-iKio + ^y^)^^^ 

since ^^ . -^-7 ^0 n -i has some small positive value while n is finite. 

0.01 (n — If (n— 2) 3rn * | r 

(37) Hence 

for 235 x . For n = <x> , z nl = — 13.5 approximately, the last figure being deter- 
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rained by the recursion formula. Computing /3 71 = 1224.4 for ^001 = — 13.5 by 
Horner's method, 

(38) ^I?»,I<0V^1?) X ( 1+ (^V') X - •*( 1 + ¥)X 122 «' 
For an infinite product, 

(40) IL,I<^ ' 2W ^ ^-n(n-y). 

(41) £001 = -8826.7. 

The first minimum is therefore a finite number. 

Second Minimum. 
The second minimum first appears in ft t = 0, and by Table IV, 

(42) iftsKl^sl, 
and both are negative. 

By equation (31) and Table IV, 

(43) ?. = (l - ^f ) &» - 0.01.25.16 ^' f ° r * = + 9 - 5j 

(44) ft, = (l - -M) As - Q01 g 6 25 & 3 , for g 78 = + 9.5; 

(45) ft 3 = (l - M) /3 73 - 001 ^ 9 36 (3 63 , for i 83 = + 9.4. 

By computing the coefficients of (43), (44), (45), it is found that (3 n3 con- 
tinues negative and of decreasing numerical value. Since the second term of 
the second member is decreasing rapidly, on account of the factors (n — l) 2 

(n — 2) 2 in the denominator, toward the limit zero, and 1 — -, — . a more 

slowly increases toward the limit unity, it is seen that ]3 m equals /3 ra multiplied 

by a positive proper fraction, plus a much smaller number. Thus (3 m , (3 m , , 

i&oo 3 remain negative and 1 ffl | > | /3 93 1 >|?oob|- 

«1™ /^«,3= /^n-1,3, and ft^ is a small negative number. By repeated use 
of (31), 

(46) £00 3 = - 27. See Fig. 4. 
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Third and Subsequent Minima. 

The third and subsequent minima can be treated in 'the same manner. We 
observe that the first minimum in any set or range may be large as compared 
with the minima in the preceding sets. For example, /3 6B , for z 6B = + 29.8, is 
— 497.6. This irregularity is due to the fact that the first minimum in a set lies 
between the last two roots of an equation of an even degree; e. g., (3 6 = 0. For 
b < 1, as in the case under consideration, the last root does not, in general, lie 
between (n — l) 2 and n 2 until b{n — 1)> 1.* 

In such a case z may be greater than (n — 1 ) 2 , as in % . Hence 

(47) ft» = (l - (^y)^« - o.oi (n— lf(n-2f & 5 haS irre g ularities - % 
computation, 

/ 29 8\ 

(48) /?„= (1 —J (-1673) -£(+3275.3), for i 6B = + 29.8. 

(49) /? 6B = (l — 1.192) (—1673) — J(+ 3275.3) = — 497.6 approximately. 

The first parenthesis changes sign, since % >(ra — l) 2 , and /? 85 therefore 
becomes a large positive number. However, the coefficient of the next term is 
comparatively large, so that the last term controls the sign of /3 65 . 

However, 

(50) & B = (l - ^) & 5 -lp a = - 152.1, for z % = 25.7 . 

Here and subsequently in /3 8B , . • ■ • , (3^ B the parenthesis is positive, slowly 
approaching the limit unity, and both terms are negative. As the last coefficient 
decreases rapidly toward the limit zero, the values of /3 8B , . . . ., (3^5 decrease 
numerically with increasing n and are always negative. Hence, as in the second 
minima, 

(51) n = °o Pn& = Pn-1, 5 

and /?qo B is a very small negative number, and not zero. 

The same argument applies to subsequent sets of minima values. 

* Heine's Kugelfunktionen, I, 407. 
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Maxima Values of (3 x . 

The first maximum of each set comes from an equation of odd degree and 
may be large, as in (3®, since z may be larger than (n — l) 2 for the last arch of 
the curve until b(n — 2) > 1 and thus make both terms of (31) positive. The 
first set of maxima values shows an increase in /3 42 on account of the small value 
of n and the factor b % [= 0.01] in the denominator of the second term. 

(52) V» = ( 1 ~ 1 ir)P«- Q.oi.4.1 P»> for 4=+H-7, 

= (— 1.9) (— 38) — 25 (— 5.85) = + 72.2 4- 146.25 = + 218.45. 

(53) ft, = (l - ^*) 0. - 0>()1 * 9 4 ft, , for z & = + 4.4 , 

= (+ 0.51) (+ 90) — 2.8 (— 87) = + 4.59 + 243.6 = + 289.5. 

(54) ff, = (l-M)^- -_l_ fl, > f OT ^ = +0.7 > 

= (+ 0.956) (4- 220) — 0.694 (— 60) = + 251.9 4- 41.6 = + 251.9 . 

The subsequent maxima (3 62 , . . • • , (3^2 will decrease and remain positive, 

1 z 

since - - -^— -— ^ decreases rapidly and 1 — -. -^ remains 

0.01 (n — l) 2 (n — 2) 2 ^ J (n — l) 2 

approximately 4- 1. 

Hence „ L i^ fi n% is a small positive number. From Table IV and equation 
(31), the approximate limit is (3^ 2 = 140. 

By the same argument, the subsequent maxima may be shown to be finite 
and not zero. 

Maxima and Minima Values of a^ . 

Max. a n = b n (n — l)! 2 . Max. /?„, (9) 

Min. a n = b n (n — 1) / 2 . Min. (3 n . 
Since the maxima and minima values of (3^, for 6 = 0.1, b = 10 and 
i = 100, are finite and not zero, 

Max. a ^ = oo . Max. (3^ = oo , 

Min. a^ = oo . Min. /^ = oo . 

Considering the factor b n in these results, it is evident that the curves 

representing a n = will be more nearly perpendicular to the z-axis for b = 10 

and 6 = 100 than for b = 0.1 . This conclusion is confirmed by Tables I, II, III 

and Fig. 5. 
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Considering the variations in a ni in these tables for z ni to four decimal places, 
it will be noticed that, in general, a ni becomes greater as n increases and the pitch 
of the curves are greater at the nullpoints. Moreover, a ni is greater in the 
vicinity of the nullpoint on the side nearest the maximum or minimum point of 
the arch, due allowance being made for the fact that the next figure of z ni in the 
table may be very large or very small in the values compared. It should also 
be noted that a point of inflexion exists in every arch, near the nullpoint which 
is more remote from the maximum or minimum point of the arch. 

Convergence op the Series 

i °° 

(55) E(q>) = , a + 2 a n cos 2n<p. 

^ n=l 

I. Argument b = 0.1 . 

The computation of a lf a z , . . . ., a n for a finite n forms a basis for more 
general conclusions. The method employed in determining the roots of a x = 
and the difficulties encountered in these computations will be illustrated by 
finding the fifth decimal figure of 2 W 3 , for argument b = 0. 1 . Having already 
z oo3— 5.5813, we compute a 6 and a 7 from (6) for 3 = 5.58138 by Horner's 
method, using the remainder to twenty-one decimal places for the determination 
of a e and a 7 . Then by the recursion formula we compute a 8 , a 9 , 

(56) a n+2 = b((n + iy~z)a n + 1 — a n . 

(57) For z = 5.58138, a 6 = + 0.22846013 

a 7 = + 0.05475009 
a 8 = + 0.00905173 
a 9 = — 0.00187210 

Evidently 5.58138 is too large for z^^ since a 9 is negative and, by (56), 
a w> a u> ■ ■ ■ • > a oo would all be negative and would increase numerically with n, 
since the coefficient of a n + 1 would be greater than 2 and would increase in- 
definitely. Hence a^z would not be zero. 

(58) For z — 5.58137, a 6 = + 0.228465748 

a T = + 0.054764586 
a 8 = 4- 0.009314650 
a 9 = — 0.000442830 

This value of z is still too great, but a 9 is much smaller numerically. 
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(59) For z = 5.58135, a e = + 0.228476975 

a 7 = + 0.054793576 
a 8 = + 0.009429440 
a 9 = + 0.000197600 
am = — 0.007939170 

Here a 9 is positive but a 10 is negative, and z = 5.58135 locates a point in the 
z-axis between the nullpoints of a 9 = and a 10 = 0. 
Computing a 1} a 2 , . . . . , a 7 by Horner's method, 

(60) For z = 5.58134, ^= — 0.279067000 

a. z = — 0.872149970 
a 3 = + 0.416988556 
a t = + 1.014702419 
a 5 = + 0.640200395 
a 6 = + 0.228482548 
a, = + 0.054808073 
a 8 = + 0.009687228 
a 9 = + 0.001783414 
a 10 = +• 0.003763039 

It is evident that 5.58134 is the value of z^ 3 , correct to five decimal places; 
since «io< a iiOi2 ••••<«<» by (56), and all the curves a 10 = to a ao = are 
still above the z-axis for z = 5.58134, but below the z-axis for z = 5.58135. 

Substituting (60) in (55), 

(61) E($) = 0.5 — 0.2790 cos 2$ — 0.8721 cos 4$ + 0.4169 cos 6$ 

+ 1.0147 cos 8$ + 0.6402 cos 10$ + 0.2284 cos 12$ 
+ 0.0547 cos 14$ + 0.0096 cos 16$ + 0.0017 cos 18$ 
+ 0.0037 cos 20$ + 

The accuracy of these coefficients is proved by the recursion formula (56). 

Convergence of the Series. 

This series is computed for a root z^g, accurate to the fifth decimal place. 
By a careful consideration of the last remainders in Horner's process used in 
computing a lf a % , . . . ., a 6 it is evident that the first three decimal figures in 
these coefficients will never change, if an infinite number of figures of z^ 3 are 
computed and substituted in E($>), since an increase of a whole unit in the sixth 
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figure of z does not cause a change in the first three decimal places of a lf a 2 , 

• ■ • • > Re- 
calling the finite sum of the first four terms F, we write, 

(62) \F($)K\F± (1.0147 + 0.6402 + 0.2284 + 0.0548 + 0.0096 + )|, 

substituting maximum values for cos 8 q>, cos 10$, . •••, and using the +sign 
when F and the following series have like signs, and the — sign when they have 
opposite signs. 

We must now define a root of a^ = more carefully. 

(63) Definition of a Boot of a^ = 0. 

From a certain a n onward indefinitel}', for an exact root of a^ = 0, a n + 1 is 
less than a n and has the same sign to n = oo . 

Otherwise one of the following relations must exist : 

(64) 1) a n + 1 ^>a n , with the same sign. 

(65) 2) a n+1 and a n differ in sign. 

(66) 3) a n + 1 <^a n for one or more terms and then a n + 1 ^>a n . 

Supposition (64) can not be true for z CX)1 a root of a^ = 0, since a, (+3 , 
a n+3> • ■ ■ • , «<» would form an increasing series, as is shown by the formula 

(67) «„ + 2 = h((n+ If — z)a n + 1 — «„• 

When b ((n + l) 3 — z) becomes greater than 2 with increasing n and 
a ll + 1 > a n , a„ + 3 must be greater than a n + 1 . Hence ^1™ a n would not be zero. 

The second supposition (65) is false for a root of a^ = 0, because (67) shows 
that, when a n + 1 and a n have opposite signs and the coefficient of a n+i becomes 
and remains positive with increasing n, a n+z must, under condition (65), have a 
larger absolute value than a n + 1 . Hence ^i™ a n is not zero. 

The third supposition (66) reduces to (64) or to (65) and therefore can not 
be true. 

Hence (63) defines a root of a^ = 0. 

This definition gives the following rule : 

(68) Rule for Computing Roots o/a co = 0. 

Find the successive figures of positive roots of a w = as great as possible, 
so that a 6 , a 7 , a 8 , . . . ., a^ shall have the same signs. 
20 
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For negative roots, each succeeding figure of z^ 1 is one less than the least 
figure that gives a permanence in the signs of a 6 , a 7 , a 8 , . . . . , a^ . 
For a root of a^ = 0, as above defined, (67) gives 

« (l+ 3 = & ((« + I) 2 — z) «n + l — ««<«« + ], 

(69) a n + ! < g-^-_ i) »_^._._ . 

From a certain n onward to n = oo , 

(70) «« + i<i«»- 
Therefore the last series in (62) may be written 

(71) 1.0147 + 0.6402+0.2284+ <1.0147+0.6402 (1 +i 4- ± + i + ....), 

which is clearly a converging series. 

Hence from (62) and (71), when z^ s is an exact root of a^ = 0, 



oo 



(72) E{ty) — s « + 2«„ cos 2n<p is finite and the series is convergent. 



»=i 



II. Argument o = 10. 

oo 

(73) E((p) = I a + 2 a„ cos 2w<?>. 



The computation of the values of a n with suificient accuracy to show that 
a n converges toward the limit zero and that the series representing E ($) is con- 
vergent when n increases without limit, involves difficulties due to the rapidity 
with which the curves a n = approach perpendicularity to the z-axis at the 
nullpoints. 

For example, we compute the third root of a 7 = by Horner's method, 
carrying the remainders to eighteen decimal places, and find the values of a 6 
and a, for this root as we obtain the successive approximations. These compu- 
tations are continued until the corresponding values of a m and a 73 give a negative 
value for a 83 by the recursion formula. 

We find z ra = 4.00133 +, and observe from the recursion formula that 

a B = 10 (49 — 4.00133) a 7 — a = 450 a 7 — a 6 , approximately. 
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Evidently a s will not become negative until 450 a 7 is less than a e 
these values of a 6 and a 7 are found, z^j can be found as follows: 



(74) 
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Z Y3 


a 73 


«63 


4 




3232143580.000000 
it 


10100579.000000 


a 


u 


1 


814029735.729124 


2543960.28 2771 


3 


88386514.036511 


276223.188814 


3 


15816887.536327 


49430.563659 


6 


1302846.520815 


4076.625965 


5 


93341.361791 


291.709237 


3 


20771.043732 


64.914147 


8 


1418.958754 


4.435519 


5 


208.453440 


0.655601 


8 


15.932590 


0.050815 


6 


1.418526 


0.005456 


5 


0.209021 


0.001676 


8 


0.015500 


0.001071 


6 


0.000956 


0.001026 


3 


0.000231 


0.001023 


9 


000013 


0.001023 


5 


0.000001 


0.001023 



Substituting the last results in (67), using z 18 to eighteen decimal places, 

(75) a 8 = 10 (49 — 4.001336538586586395) . 0.0000012 — 0.0010230 

= _ 0.000490. 

Since a 8 is negative, z is too large by (63). 
Taking z to seventeen decimal places, 

(76) a 8 = 10 (49 — 4.00133653858658639). 0.0000133 — 0.0010231 

= + 0.0049617. 

Since a s is here positive and in the former case negative, z co8 lies between 
the two values taken. 

2o0 3 = 4.00133653858658639 +. 

Subsequent figures may be found as in (57) to (60). 
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Hence, 

5 

(77) E($) = 0.5+ 2 a„cos2w0+O.OOlO231 cos 120+0.0000013 cos 140+ 



n=i 



The trial divisor in Horner's process contained twenty figures when the 
contract method was used in determining the last fifteen figures of z. To secure 
values for a 1 to a s that will satisfy the recursion formula, twenty-one decimal 
places should be used in all remainders and the contract method should be em- 
ployed later in the work. 

It will be noted that even with seventeen decimal places of z M 3 accurately 
computed, a^ is greater than a 73 . It is evidently practically impossible to find 
a value for 2^,3 so near the exact root that a 73 , a 83 , ...., a^ 8 shall form 
a rapidly decreasing series. Heine* remarks that when such a value is found, 

(78) «»+i< 



6[(n + l) 8 — z] — 1' 

and therefore with increasing n, a n decreases rapidly to the limit zero Practically 
each succeeding figure of z cai must be found as in (57) to (60). 

For argument b = 100, similar treatment will give like conclusions. As 
the pitch of the curves is greater at the nullpoints, more decimal places will be 
required and greater difficulties in computation will be encountered. 

The Best Values of b and z. 
To find a finite number of terms of 

(79) ^(4>) = i a o + a i cos 20 + a 2 cos 40 + . . . . 

that will give a good approximate solution of the equation 

( 8 °) Cl ^ + (1 C0S 2 * + 4Z ) E{<p) = °> 

it is evident from the preceding discussion that z M e must be computed to four or 
more decimal places when b = 0.1 and to many more places when & is a larger 
number. It will be seen that small values of b are desirable ; since the curves 
representing a n = have a steeper pitch for larger values of b, and consequently 
a larger number of decimal figures of z^ t must be computed in the latter case 
to obtain a n and a n + 1 small enough to satisfy the necessary condition given in 

* Kngelfmiktionen, I, 411. 
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equation (4). Fairly good values of the coefficients of our series for b = 0.1 
have been found in computing z^ j to z^ 7 . 

5 

(81) 2001 = — 16.9015. ^ = 0.5 + 2 a n cos 2n$ + 0.004242 cos 12$ 

» = i 

+ 0.003259 cos 14$ + 

5 

z CX)2 =— 5.0524. E 2 = 0.5 + 2 a n cos 2n$ — 0.055101 cos 12$ 

— 0.013109 cos 14$— 0.00 246 9 cos 16$— 0.001497 cos 18$— 

5 

ZaoS =4- 5.5813. # 3 = 0.5 4- 2 a n cos 2rc$ + 0.228301 cos 12$ 

» = i 

4- 0.054866 cos 14$ +0.009922 cos 16$ + 0.003074 cos 18$+ 

5 

Zoo i = + 14.5616. E t = 0.5 + 2 a n cos 2n$ — 0.434188 cos 12$ 

n = 1 

— 0.131937 cos 14$ — 0.030213 cos 16$— 

5 

Zoo 5= + 20.4705. # B =:0.5+ 2 a n cos 2«$ + 0.627187 cos 12$ 

re=l 

+ 0.233685 cos 14$ + 0.039501 cos 16$ + 0.005334 cos 18$ 

+ 0.002985 cos 20$+ 

5 

3^6= + 27.1931. #0 = 0.5+ 2 a n cos 2»$ — 5.360337 cos 12$ 

— 2.830106 cos 14$ — 0.811261 cos 16$ — 0.155902 cos 18$ 

— 0.027635 cos 20$— 

6 

2^,= + 37.3476. #7 = 0.5+ 2 a n cos 2?i$ + 69.726061 cos 12$ 

n — l 

+ 51.778903 cos 14$ + 28.841160 cos 16$ + 6.882976 cos 18$ 

+ 1.142740 cos 20$ + 0.166280 cos 22$ + 

It will be noticed in these equations, that, in general, the smaller the roots 
are algebraically the better are the coefficients obtained for an approximate 
solution in the form of a finite number of terms of the infinite series. From 
Fig. 5, we should expect the best approximations for # 5 and # 4 ; and this is 
doubtless in general true, since the maxima and minima are here less numerically 
and the slopes of the curves at the nullpoints are not so great, but the possible 
difference in the fifth and subsequent decimal places of %&}, . ■ •-, z s! make 
it impossible to determine this fact definitely without further calculations. 

To obtain, for b = 10,. as good approximate results as the above, the values 
of 2ooi, ••••, Zoo, must be computed to eighteen or more decimal places and, 
for b = 100, to many more places. 



152 Butts: The Elliptic Cylinder Function of Glass K. 

As a general conclusion, it is evident that small values of b and z are the 
best values. 

General Proof of Convergence. 

Granting that a^ =zf(b, z) = is an equation of an infinite degree in the 
general form of equations (5) and that the coefficients of the series in (2) satisfy 
the recursion formula (3) for exact roots of a^ = 0, we have shown that these 
roots can be computed to any desired degree of accuracy by rule (68). We must 

00 

now show that the series in the solution JE(<p) = ia + 2 a n cos 2n<p is con- 

« = i 

vergent for values of a n computed for 2 W { , any exact root of a^ = 0. 
Beginning with a certain a m , definition (63) gives 

«m + l<i«m- (70) 

00 

(82) E(ty) = \ a + 2 a n cos 2n<p 

71 = 1 

— (1 a o + a i cos 2$ + • • • • + a m -\ cos 2{m — 1) ty) 
+ (a m cos 2m$> + ■••• + a oo cos °° #)• 
If a m is the first term of the decreasing series that characterizes a root of 
a oo = in definition (63) and satisfies relation (70), 

I E W) |< I (i «o + «i cob 2$ + + a m _j cos 2 (m — 1 ) $) 

± (a m + a ro + i + + a a0 )\, 

substituting maximum values for cos 2m<p, . . . . , and using the plus sign when 
the two parentheses are alike in sign, and the minus sign when they are of 
opposite sign. 

|«m + «m+l + ----+«oo l<l«m(l + i + i + ••••) 1 = 1 2a m|- (70) 

The quantities a 1; a 2 , . . . ■ , a m are finite, since a n is a rational, integral 
function of z with finite coefficients; hence a 1} a 2 , ...., a m must be finite for 
finite values of z. 

Hence, 

I E{$) l< ! (Finite Sum ± 2a m ) \ , and 



E(<p) = i a + 2 a n cos 2n$ is a convergent series and a solution of 
« + (|cos2*+4 z )i?=0. 
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TABLE I. 
b = 0.1 



TABLE II. 
6 = 10 



TABLE III. 
b = 100 



Z n 


+ 0.0000 


A n 


+ 0.00000 


Z n 


+ 0.0000 


A n 


+ 0.00000 
xlOi 


z n 


+ 0.0000 


A n 


+ 0.00000 

XlO« 


Z.1 


-13.6509 


A 2i 


— 0.00000 


Z 21 


- 0.0196 


A 21 


- 0.00002 


Z 21 


- 0.0001 


21 


- 0.00004 


z 22 


+ 14.6509 


A 22 


-0.00000 


Z 22 


+ 1.0196 


A 22 


- 0.00002 

X102 


Z 22 


+ 1.0001 


A 22 


- 0.00004 
xlO« 


Z 31 


-16.2479 


A 31 


— 0.00002 


Z S1 


- 0.0196 


A 31 


- 0.00013 


Z 31 


- 0.0001 


31 


- 0.00019 


Z 32 


+ 3.6470 


A 32 


—0.00000 


Z 32 


+ 1.0163 


A 32 


- 0.00006 


Z 32 


+ 1.0001 


32 


- 0.00009 


Z 33 


+ 18.6009 


A 33 


+ 0.00002 


Z 33 


+ 4.0033 


A 33 


+ 0.00071 
xlO» 


Z 33 


+ 4.0000 


A 33 


+ 0.00020 
xl06 


z« 


-16.8064 


A 41 


— 0.00000 


Z 41 


- 0.0196 


A 41 


- 0.00112 


Z 41 


- 0.0001 


A 41 


— 0.00179 


z « 


- 2.8853 


A 42 


+ 0.00002 


Z 42 


+ 1.0163 


A 42 


- 0.00032 


Z 42 


+ 1.0001 


A 42 


- 0.00079 


z 43 


+ 12.3932 


A 43 


+ 0.00000 


Z 43 


+ 4.0013 


A 43 


+ 0.00110 


Z 43 


+ 4.0000 


43 


+ 0.00040 


z« 


+ 21.2985 


A 44 


-0.00001 


Z 44 


+ 9.0019 


A « 


- 0.01395 
xlO* 


Z 44 


+ 9.0000 


A 44 


- 0.00359 
x 10 8 


Z 5. 


-16.8934 


A 51 


—0.00015 


Z 51 


- 0.0196 


A 51 


- 0.01875 


Z 51 


- 0.0001 


i 51 

■a-52 


— 0.02879 


z 52 


— 4.6582 


A 52 


+ 0.00003 


Z 52 


+ 1.0163 


A 52 


- 0.00476 


Z 52 


+ 1.0001 


— 0.01199 


Z 53 


+ 7.7970 


A 53 


+ 0.00001 


Z 53 


+ 4.0013 


A 53 


+ 0.01315 


^<i3 


+ 4.0000 


t 


+ 0.00480 


z 54 


+ 18.0760 


A 54 


-0.00002 


Z54 


+ 9.0005 


A 54 


- 0.09031 


z 


+ 9.0000 


- 0.00719 


Z 55 


+ 25.6784 


A 55 


+ 0.00034 


Z 55 


+ 16.0014 


A 55 


+ 0.28200 
x 105 


< 


+ 16.0000 


A 55 


+ 0.14400 
xlO 10 


Z 6, 


-16.9011 


A 6! 


-0.00013 


Z 61 


— 0.0196 


A 6I 


- 0.52699 


Z 61 


- 0.0001 


A 61 

A 62 
A 63 
A 64 

A 5 

66 


— 0.71997 


Z 62 


- 5.0133 


A 62 


+ 0.00008 


Z 62 


+ 1.0163 


A 62 


— 0.08015 


^62 


+ 1.0001 


— 0.28791 


Z 63 


+ 5.9979 


A 63 


+ 0.00004 


Z 63 


+ 4.0013 


A 63 


+ 0.27626 


^63 


+ 4.0000 


+ 0.10080 


Z 64 


+ 15.7525 


A 64 


-0.00002 


Z 6 4 


+ 9.0005 


A 64 


— 1.44069 


Z 64 


+ 9.0000 


- 0.11520 


Z 6 5 


+ 31.9884 


A 65 


+ 0.00003 


Z 6 5 


+ 16.0003 


A 66 


+ 1.59149 


Z 65 


+ 16.0000 


+ 0.28798 


Z 66 


+ 33.1754 


A 66 


-0.00074 


Z 66 


+ 35.0011 


A 


—10.02866 


7 


+ 35.0000 


— 8.83001 
















X106 








x 10 12 


Z 71 


— 16.9015 


A 71 


— 0.00109 


Z 71 


- 0.0196 


A 71 


— 16.90346 


Z 71 


- 0.0001 


A 71 


-25.91922 


Z 72 


- 5.0500 


A ,2 


+ 0.00008 


Z 72 


+ 1.0163 


A ,2 


- 3.86381 


Z 2 


+ 1.0001 


A 

A 73 

A 74 

I" 
A 76 

A 77 


-10.07714 


Z 73 


+ 5.6193 


A 73 


+ 0.00000 


Z 73 


+ 4.0013 


A 73 


+ 3.33886 


Z 73 


+ 4.0000 


+ 3.22567 


Z 74 


+ 14.7525 


A 74 


-0.00004 


Z 74 


+ 9.Q005 


74 


- 7.36793 


Z 74 


+ 9.0000 


- 3.11034 


Z 75 


+ 30.7688 


A 75 


+ 0.00006 


z„. 


+ 16.0003 


A 75 


+ 31.69049 


Z 75 


+ 16.0000 


+ 5.76028 


Z 76 


+ 28.6580 


A 7 6 


— 0.00021 


z ™ 


+ 25.0002 


A 76 


— 196.55479 


Z 70 


+ 25.0000 


— 20.16004 


Z 77 


+ 43.1526 


A 77 


+ 0.00348 


Z 77 


+ 36.0009 


A 7, 


+ 1084.83159 


Z 77 


+ 36.0000 


+ 1088.63940 
















xlO' 








x 10" 


Z «>1 


— 16.9015 






Z oo 1 


— 0.0196 






Z aol 


- 0.0001 






Z 0O2 


— 5.0524 






Z 0O 2 


+ 1.0163 






Z Q0 2 


+ 1.0000 






Z 00 3 


+ 5.5813 






Z OD 3 


+ 4.0013 






Z X 3 


+ 4.0000 






Z 0O 4 


+ 14.5616 






Z CO 4 


+ 9.0005 






Z 00 4 


+ 9.0000 






Z CO 5 


+ 20.4705 






Z CO 5 


+ 16.0003 






Z 00 5 


+ 16.0000 






Z 0O6 


+ 27.1931 






Z 0O 6 


+ 25.0002 






Z Q0 fi 


+ 25.0000 






Z 0O 7 


+ 37.3476 






Z 00 7 


+ 86.0001 






7 
^CO 7 


+ 36.0000 
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TABLE IV. 



TABLE V. 





6 = 0.1 




+ p = W 


ax. 




6= 10 


+ P = 


Max. 








— ft = Min. 








~P = 


Min. 


a n - 


= 0, /? n =0 


z 


P 


a 


a n = 


= o, p n = 


2 


J8 


d 


Z 21 


-13.6509 








Z 2, 


- 0.0196 












+ 0.5 


- 100.1 


- 1.00100 






+ 0.5 


— 0.13 


- 0.00130x10* 


7 

' J 22 


+ 14.6509 








Z 22 


+ 1.0196 








Z 31 


-16.2479 








Z 3. 


— 0.0196 












- 8.4 


- 292.6 


— 1.17040 






+ 0.4 


- 0.11 


- 0.00048x106 


Z 32 


+ 2.6470 








Z 32 


+ 1.0163 












+ 11.7 


+ 218.2 


+ 0.87280 






+ 2.8 


+ 0.75 


+ 0.00300x106 


Z 33 


+ 18.6009 








z 


+ 4.0033 








Z 41 


-16.8064 








z « 


- 0.0196 












-11.4 


— 497.4 


— 1.79071 






+ 0.4 


- 0.10 


- 0.00036 xlO 8 


Z 42 


- 2.8853 








Z 42 


+ 1.0163 












+ 4.4 


+ 289.6 


+ 1.04256 






+ 2.7 


+ 0.52 


+ 0. 00187 x10 s 


Z 43 


+ 12.3932 








Z 43 


+ 4.0013 












+ 17.8 


- 188.0 


— 0.67680 






+ 7.7 


— 3.43 


- 0.01235 xlO 8 


Z 44 


+ 21.2985 








Z 44 


+ 9.0019 








Z 5! 


— 16.8934 








Z 51 


- 0.0196 










-12.6 


- 708.8 


- 4.08311 






+ 0.4 


— 0.10 


— 0.00062 xlO 10 


Z 52 


- 4.6582 








Z 52 


+ 1.0163 










+ 0.7 


+ 252.4 


+ 1.45433 






+ 2.7 


+ 0.45 


+ 0.00261xl0"> 


Z 53 


+ 7.7970 








Z 53 


+ 4.0013 










+ 12.9 


- 156.7 


- 0.81259 






+ 7.0 


- 1.96 


— 0.01134 xlO 1 " 


Z 54 


+ 18.0760 








Z 54 


+ 9.0005 










+ 22.7 


+ 192.9 


+ 1.11110 






+ 13.7 


+ 15.84 


+ 0.09129 xlO 10 


Z 55 


+ 25.6784 








Z 55 


+ 16.0014 








Z 6! 


-16.9010 








Z 61 


- 0.0196 










-13.0 


— 1137.2 


— 16.37613 






+ 0.4 


- 0.10 


- 0.00153 xlO 12 


Z 62 


- 5.0133 








Z 62 


+ 1.0163 










- 0.5 


+ 205.5 


1- 2.96017 






+ 2.6 


+ 0.38 


+ 0.00558 xlO 12 


Z 63 


+ 5.9979 








Z 63 


+ 4.0013 










+ 9.5 


- 86.1 


- 1.24016 






+ 6.9 


— 1.41 


- 0.02041 xlO 12 


Z 64 


+ 15.7525 








Z 64 


+ 9.0005 










+ 19.2 


+ 49.1 


+ 0,70733 






+ 13.4 


+ 7.20 


+ 0. 10364 xlO^ 2 


Z 65 


+ 21.9884 








Z m 


+ 16.0003 










+ 29.8 


- 497.6 


- 7.16590 






+ 22.3 


-69.33 


- 0.99847 xlO' 2 


Z 66 


+ 33.1754 








Z 66 


+ 25.0011 








Z ,l 


-16.9015 








Z 71 


- 0.0196 










-13.3 


-1227.4 


-63.63132 






+ 0.4 


- 0.10 


- 0.00547x10" 


Z 72 


— 5.0500 








Z 72 


+ 1.0163 










- 1.0 


+ 165.9 


+ 8.60357 






+ 2.6 


+ 0.35 


+ 0.01850x10" 


Z 73 


+ 5.6193 








Z 73 


+ 4.0013 










+ 9.1 


- 47.0 


- 2.44122 






+ 6.9 


- 1.14 


- 0.05941x10" 


Z ,4 


+ 14.7525 








Z U 


+ 9.0005 










+ 17.8 


+ 23.2 


+ 1.20588 






+ 13.2 


+ 4.51 


+ 0.23393 x 10" 


z_. 


+ 20.7688 








Z 75 


+ 16.0003 








,0 




+ 25.7 


- 152.1 


- 7.87878 






+ 21.9 


-26.29 


— 1.36291x10" 


Z 76 


+ 28.6580 








Z 76 


+ 25.0002 












+ 39.5 


+1431.7 


+ 74.22189 






+ 34.0 


+ 1551.99 


+ 80.45630x10" 


z n 


+ 43.1526 








Z»7 


+ 36.0009 









a„ = 6» (»- l)/ 2 A, 



a = 6" (» - 1) / 2 3 
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